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ABSTRACT 

We investigate the cross-correlation between the cosmic infrared background (CIB) 
and cosmic microwave background (CMB) anisotropics due to the integrated Sachs- 
Wolfe (ISW) effect. We first describe the CIB anisotropics using a linearly biased power 
spectrum, valid on the angular scales of interest. From this, we derive the theoretical 
angular power spectrum of the CMB-CIB cross-correlation for different instruments 
and frequencies. Our cross-spectra show similarities in shape with usual CMB/galaxies 
cross-correlations. We discuss the detectability of the ISW signal by performing a 
signal-to-noise (SNR) analysis with our predicted spectra. Our results show that : 
(i) in the ideal case of noiseless, full-sky maps, the significances obtained range from 
6 to 7a depending on the frequency, with a maximum at 353 GHz (ii) in realistic 
cases which account for the presence of noise including astrophysical contaminents, 
the results depend strongly on the major contribution to the noise term. They span 
from 2 to 5a, the most favorable frequency for detection being 545 GHz. We also find 
that the joint use of all available frequencies in the cross-correlation does not improve 
significantly the total SNR, due to the high level of correlation of the CIB maps at 
different frequencies. 



1 INTRODUCTION 



The discovery of the acceleration of the expansion of the 
Universe, made through supernovas observations (Riess 



et al. 1998 Perlmutter et al. 1999 1 at the end of the last 



century, has since led to many theories aimed at explaining 
its origin. These theories have been regrouped under the 
term "Dark Energy" (DE), designating a new and unknown 
component of our Universe which theoretically amounts to 
70% of its total energy budget. Among the many solutions 
proposed to account for this intriguing phenomenon, one 
of the leading contenders is the so-called "cosmological 
constant", an idea first introduced by Einstein in his 
original theory of general relativity to achieve a stationary 
universe, but which he discarded after the discovery of the 
Hubble redshift. This cosmological constant is assimilated 
to an intrinsic energy density of the vacuum, and therefore 
is constant in time and space ; it also has an equation of 
state w = p/ p equal to —1, both on theoretical grounds and 
because no confirmed deviations from w — —1 have been 
detected so far. Despite its simplicity, it does reproduce 
most of the current observations while being (one of) the 
most "economical" solution, but it is nevertheless plagued 
by a few serious theoretical problems (e.g. [Padmanabhanj 
2003{|Bass|20Tl| ). 



Apart from these theoretical issues, the accelerated 
expansion of the Universe still needs to be tested further 
in the framework of the A-CDM model by independent 
measurements from cosmological observations. Over the 
last decade, other possible probes have been proposed such 
as the study of baryon acoustic oscillations (Eisenstein, Hu 
fc Tegmarkl|1998| [Eisenstein et al . 2005; Bassett & Hlozek 



2010 



and references therein) which provide a "standard 
in cosmology, and are heavily influenced by the 
and so, by the dark energy. 



ruler' 

energy content of the Universe 



In this article, we focus on an alternative probe of 
the Dark Energy, namely the study of the integrated 
Sachs-Wolfe (ISW) effect ; the "or iginal" SW effect first 
introduced at the end of the 60s (Sachs & Wolfe 19671 



describes the imprint on the cosmic microwave background 
(CMB) of anisotropies caused by gravitational redshift 
occurring at the surface of last scattering. Its "integrated" 
counterpart is similar in that it also has a gravitational 
origin and contributes to the CMB secondary anisotropies, 
but it only occurs in a Universe not dominated by matter. 
Indeed, the ISW effect is caused by the large-scale structures 
of the Universe, whose gravitional potentials are slowly 
decaying - instead of being constant in a matter-dominated 
regime - and therefore giving a net gain (in case of a 
potential well) or loss (hill) of energy to the CMB photons 
that travel across them. 



© 0000 RAS 



2 Ilic, Douspis, Langer, Penin & Lagache 



This effect shows in the power spectrum of the CMB 
temperature anisotropies at large angular scales (low £) but 
the cosmic variance at those very multipoles together with 
the relatively small amplitude of the ISW effect make its 
direct detection very challenging, if not impossible, when 
using only the CMB itself. To circumvent this limitation, 
cosmologists have devised a way to exploit the link between 
this imprint on the CMB and the large-scale structures 
causing it, by simply cross-correlating the CMB with matter 
density maps (actually galaxy maps in practice) and then 
comparing the results to a null hypothesis and to what is 
expected from theory. 

During the last decade or so, a growing interest has 
arisen in this field thanks to the development of large 
galaxy surveys : SPSS ([Abazajian et al.|2009[ for the latest 
releas e), NVSS ( [Condon et al.|1998[ ), 2MASS ( [Jarrett et al. 
2000), etc. They allow cosmologists to cross-correlate the 



CMB (as seen by WMAP) to proxies of the matter density 
field as seen at many wavelengths : X-rays ( Boughn & Crit- 
tenden||2004| ), optical ([Granett, Neyrinck fc Szapudi||2009f7 
near-infrared ([Rassat et al.|2007 l or radio (Pietrobon, Balbi 
& Marinucci 2006). However, this method has yet to pro- 



duce a definitive and conclusive detection of the ISW effect. 



with significances so far ranging from negligible ( Sawangwit 
et al.||2010[ ) to 4.5o- ( [Giannantonio et al.|[2008| throughout 
the literature. The potential of future surveys such as 
LSST, Pan-STARRS or Euclid, have also been explored in 



terms of signal-to-noise ratio of the ISW detection ( Douspis 
et al.|2008 |. Another noteworthy approach by Taburet et al. 
(20101 considered the cross-correlation of the ISW effect 



with the thermal Sunyaev-Zel'dovich effect as both effects 
take place in the same potential wells ; this could provide 
an independent probe for the existence of Dark Energy out 
of pure CMB data. 

The originality of our work is to consider here the 
Cosmic Infrared Background (CIB), first discovered by 
[Puget et al.| ( |1996| ). This background, visible roughly from 
10 to 1000 fj,m in wavelength, arises from accumulated 
emissions from star-forming galaxies spanning a large range 
of redshifts. The earliest epoch for the production of this 
background is thought to be when star formation first 
began, and contributions to the CIB continued through the 
present epoch, including our current dark energy dominated 
era. The CIB also features anisotropies (first detected and 
discussed in Lagache fc Puget|2000 Matsuhara et al.|2000 1 
that are underlined by the galaxy density field and thus 
the matter density fiuctuations. It is therefore reasonable 
to expect that it has a positive correlation with the CMB 
through the aforementioned ISW effect. 

In this paper, we first present an analytical calculation 
of the CMB-CIB cross-correlation signal through the 
integrated Sachs-Wolfe effect. We then use it to compute 
the expected power spectrum of this correlation in different 
cases, namely with a CIB observed at several frequencies 
and with various instruments (IRAS, Herschel-SPTKE, 
and Planck-BFl) . With these results we perform a signal- 
to-noise ratio (SNR) analysis in order to quantify the 



detectability of the cross-correlation, focusing first on the 
"perfect case" scenario, i.e. a situation where both the 
CMB and CIB are full sky maps, without noise, so that 
the detection is only limited by the cosmic variance. We 
then discuss the effect of noise (including contaminating 
astrophysical components and instrumental noise) in the 
maps and its consequences on the SNR. Finally, we end with 
a few conclusions and discussions about the perspective of 
application of our predictions. 

Throughout all our calculations we assume a Euclidean 
Universe corresponding to the WMAP 7 best-fit cosmology, 
with adiabatic scalar perturbations and a nearly scale in- 
variant initial power spectrum. 



2 MODELLING THE EXPECTED SIGNAL 
2.1 CIB anisotropies 

Ever since its discovery, many efforts have been deployed 
to detect the cosmic infrared background with increasing 
precision, especially in order to study its anisotropies which 
contain a lot of information about the star and galaxy 
formation histories, including their clustering processes. 
The most recent papers on the CIB anisotropies use 
sophisticated models which compute the halo occupation 



distribution (HOD, see e.g. Peacock fc Smith||2000 Cooray 



fc Sheth 2002 1 and the Dark Matter halos properties, in 



order to predict the power spectrum of these anisotropies. 
Recently applied to the new Planck data (sec Planckl 
Collaboration|2011b I, this framework allowed us to confirm 
that the bias between infra-red galaxies and the linear 
theory matter power spectrum is not independent of scale, 
and that the halo occupation distribution is evolving with 
redshift. 

Such models are particularly useful when describing the 
small, non-linear scales of the CIB. Since we focus here on 
the ISW effect which only concerns much larger scales, we 
can use a simpler model for the CIB, similar to the descrip- 



tion made by Knox et al. (20011. The general definition of 



the CIB anisotropies at a given frequency u and in a given 
direction n can be then written as the following line-of-sight 
integral : 



STciB{n,u) = 



^0 ^„ 

dz — a{z)5j{{T]o-ri)h,u,z) (1) 



with Sj being the emissivity fluctuations of the CIB. The 
integration is made over rj, the conformal time, from some 
initial time rj;^,- before star formation began to our location 
at the coordinate origin -qo- In their work, Knox et al. hy- 
pothesized that the CIB anisotropies are direct tracers of the 
matter density fiuctuations 5 — Spm/Pm, up to a bias fac- 
tor. Therefore, the previous expression becomes an integral 
of the product between a mean far infrared (FIR) emissivity 
and the matter density fluctuation fleld : 

STciB{n,v)= / dz — a{z)bj(iy,z)j{u,z)S{{r]o-ri)h,z) 

(2) 
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Figure 1. Angular power spectra of the CIB fluctuations at four frequencies of tiie Piancfc-HFI instrument, as predicted by the Planck 
Team (blue continuous line) and by our non- biased models (dashed yellow line). For each frequency, we provide in red the linear bias 
which gives the best agreement between the two models, and plot our models taking into account this bias (solid yellow line). The data 
points correspond to measurements obtained by the Planck Team | |Planck Collaboration|2011b[ |. 



The quantity bj (v, z) is a frequency- and redshift-dependent 
matter-emissivity bias defined by : 



= bj{u,z) 5{{r]o-v)A,z) (3) 



and j(i/, z) is the mean emissivity per comoving unit volume 
at frequency f as a. function of redshift z, which is derived 
here using the empirical, parametric model of |Bethermin| 
et al. (20111. The matter density field 5 is described in our 



analysis by a linear power spectrum. While this approxi- 
mation is not accurate at small scales where non-linearities 
arise, it is perfectly valid for the scales of interest in our 
work since the first hundred multipoles {i < 100) comprise 
most of the ISW signal. Following a calculation similar to 
the one described in the next section, we can express the 
angular power spectrum of the CIB fluctuations as follows : 



^CIB/ X , 9 



dk 



(4) 



where Mi{k, u) is given in Eq. ( 10 1 and is defined below 
Eq. ([8}. 



Lastly we choose the previously mentioned linear bias[^ 
bj{i',z) to be constant in redshift : bj{i',z) — 6iin(//). To 
obtain it at each frequency, we compute the value of bun 
that gives the best agreement between our linear CIB power 



spectrum and those obtained from the Planck data ( Planck 



Collaboration 2011b I. We choose to fit the two spectra in the 
range of multipoles £ G [10, 50], where most of the ISW sig- 
nal resides. This is illustrated in Fig. [l] where we plot the bi- 
ased and non-biased CIB linear spectra from our framework 
and compare them to the ones from [Planck Collaboration! 
(2011b I at their four frequencies. Overall, the two sets of 



spectra show good agreement over the multipoles of interest 
; the spectra deviate at higher £ (starting from ~ 100) due to 
the rise of non-linearities that we did not account for in our 
linear model - namely the small-scale correlations between 
galaxies inside the same halos. The linear bias we obtain this 
way increases with the wavelength : this is coherent with the 
fact that as we go further deep into the infrared, the galaxies 



^ This bias here represents our matter-emissivity bias in Eq. |3| 
and should not be confused with the widely used galaxy-Dark 
Matter bias, though ours does contain information about how 
the emitting objects populate Dark Matter halos. 
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Figure 2. Theoretical angular cross power spectrum of the CIB-CMB correlation calculated for IRAS at 100 iim (left-hand panel), 
for Herschel SPIRE between 250-500 fim (central panel) and Planck HFI between 350-2097 fim (right panel). The linear bias, bun, is 
fixed here to 1 at all frequencies in order to compare the non-biased CIB power spectra. The vertical dashed line on each panel marks 
the upper limit of the multipoles used in our analysis : this choice comes from the absence of ISW signal (see Fig. [sjl and the rise of 
non-linearities at higher £. 



probed are more luminous at higher z. They reside in more 
massive and rarer halos, and are therefore more biased. 

2.2 Correlation with the ISW 

In the CMB anisotropies, the temperature contribution due 
to the ISW effect is an integral over the conformal time of 
the growth rate of the gravitational potentials : 



5Tisw(h) = / dri e 



-T(rj) 



(<i>-*)[(r?o-r?)n,»7] (5) 



where rir is some initial time deep in the radiation era, <I> 
and ^ are the Newtonian gauge gravitational potentials 
(with the conventions used in Kodama fc Sasaki|1984 1, T{ri) 
is the optical depth included to account for the possibility 
of late reionisation, and the dot denotes differentiation with 
respect to r;. 

We are interested in calculating the CIB-ISW cross- 
correlation function C"^"^ at a given frequency u in direct 
space : 

C"(efli,fi2,i/) = {STciB{Ai,u)5TiswiA2)}. (6) 
After a decomposition into Legendre series, we get : 



C"(0, = ^ Cr(t.) A(cos e) 



4n 



(7) 



where we do not include the monopole and dipole terms in 
the sum. Using Eqs. (pl and ([5|, we follo w a calculation 
similar to [Garriga, Pogosian fc Vachaspati] ( |2004[ ) in order 
to finally get the CMB-CIB cross-power spectrum at a fre- 
quency u : 

where comes from the primordial curvature power 

spectrum P-r. = 27r^A^/fe^. The use of this primordial 



spectrum differs from previous works on CMB-galaxies 
cross-correlation, where the present matter power spectrum 
is usually introduced instead, and is then evolved backward 
in order to find its correlation with the CMB. Conversely, 
in the Garriga et al. approach, the starting point is the 
primordial perturbations which are evolved to the present 
time. While it allows a full account of possible fluctuations 
in the dark energy in non-A models, it also avoids the 
frequently used Limber approximation, which is known to 
be somewhat inaccurate at the largest scales considered 
here. 

At this point, we need to compute the two main func- 
tions Tf"^ and Me, which are defined as : 

rf ^(fc) = r'dn e'"*"' j,{k[n - no]) {c^^i> ~ 4>), (9) 



and 



Mt{k, v) = / d-q ji {k[ri - ryo]) a{r]) biin{u) j{u,ri) S{k, rj) 
•I vo 

(10) 

where je{-) are the spherical Bessel functions, while 5, (j> and 
■(/; are the time-dependent parts of (respectively) the Dark 
Matter density contrast 5, and the two Newtonian gravita- 
tional potentials $ and 5". The two coefficients and cs<s/ 
give the relations between 5, <E> and for adiabatic initial 
conditions : 



-5 3 _* A2„ 

= ^ = "2 ' = $ = " i ^ + 5^'^ 



(11) 



where = pv/{pv + p~,), pv and p^ being respectively the 
energy densities in relativistic neutrinos and photons. 



^ The separation between time and space dependance in the 
terms <5, "J" and is allowed in our calculation since the time 
evolution of each Fourier mode only depends on the magni- 
tude k = ||k||. For exemple : <3?(k, ry) = <I>(k, r7r)0(fc, J?) ■ 
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Figure 3. Left-hand panel : Cumulated SNR as a function of imax (defined in Eg. |12| for tiie CMB-CIB cross-correlation, at our chosen 
frequencies and instruments. Right-hand panel : Total SNR with imax = 100 as a function of frequency/wavelength. 



2.3 Shape of the cross-correlation spectrum 

To compute these expressions, we adapted for our analysis 
an already modified version of CMBFAST ([Seljak fc 
Zaldarriaga]|1996[), named CROSS -CMBFAST ( |Cora,saniti," 
Giannantonio fc Melchiorri| |2005| . For a given cosmology 
and emissivity function j{v,z) (see Eqs. |2| and (|3|), our 
code calculates the Cf from Eq. Q and at the same 
time the predicted power spectrum of the GIB fluctuations 
described in Eq. Q and already illustrated in Fig. [l] It also 
gives the standard GMBFAST outputs, including the GMB 
temperature power spectrum. 

In Fig. [2] we present our predictions for the GIB-GMB 
cross-correlation, at several FIR wavelengths and for 
different instruments, namely : IRAS at 100 /xm, Herschel 
SPIRE at 250, 350 and 500 /xm and Planck HFI at 350, 
550, 850, 1380, and 2097 /xm. We note that at 350 the 
SPIRE- and Planck- predicted spectra differ slightly from 
each other, due to the difference in wavelength bandwidth 
of the two instruments. 

In a fashion similar to previous galaxy-ISW cross- 
correlations (see the references in Section [TJ, we note that 
the cross-correlation peaks around I ~ 10-30, and quickly 
vanishes at higher multipoles. Gomparing the signal at the 
different wavelengths shows that the amplitude of the cross- 
correlation signal is maximum at a wavelength ~ 250 /im. 
This is not entirely surprising, since this wavelength roughly 
corresponds to the maximum of the observed GIB spectral 
energy distribution (see Dole et aL]|2006 for reference). 



It should be also noted that these results are not exact 
at the highest Is since the non-linear counterpart to the ISW 
effect, called the Rees-Sciama effect, contributes at those 



discussion). However, in our case the linear part of the ISW 
largely dominates at the observed peak in Fig. [2] 



3 SIGNAL-TO-NOISE ANALYSIS 
3.1 Ideal case 

We now investigate the detection level of the ISW effect 
using GMB-GIB cross-correlation by performing a signal-to- 
noise ratio analysis. Using the power spectra computed in 
the previous section, we can write for each given frequency 
1/ the total signal-to-noise ratio of the ISW detection as : 



[C-(//)]2-HCCIB(l.) xCfMB 



(12) 



scales (see Schaefer, Kalovidouris fc Heisenberg^2010j for a 



where the total (or cumulative) signal-to-noise is summed 
over multipoles between 1 = 2 and ^max 5; 100 where the 
signal has its major contribution (see previous section, 
Fig.[2f. 

In this section, we first consider the ideal situation 
where the GIB and GMB maps used for cross-correlation are 
noiseless and cover the whole sky ; with these assumptions 
we obtain the highest possible signal-to-noise ratio, the only 
limitation being the cosmic variance. In Fig. |3] we present 
our prediction for the GIB-GMB cross-correlation in the 
case of a full-sky GIB map, providecj^ by the previously 
mentionned instruments and frequencies. 



^ Only the IRAS 100 ^im data is already available, and previ- 
ous works have managed to extract its CIB component on small 
patches of sky ( [Miville-Deschenes, Lagache, fc Puget|2002[ |, but 
the CIB has yet to be extracted over a large enough part to allow 
for an ISW detection. 
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With these optimistic assumptions, we obtain high 
levels of detection for the CIB-CMB correlation which reach 
~ 7(7 (for detailed results, see Table [T]). It should be men- 
tioned that these results in the ideal case are independant 
of the previously discussed linear bias in Section |2.1| even 
if it boosts the correlation signal. This can be understood 
from Eq. ( |12[ | where the linear bias can be factorized from 
each term (one for C" and a squared one for Cf^^) and 
therefore cancels out. 



CIB, therefore reducing the "usable" fraction of the sky by 
at least ~ 25%. Furthermore, the rest of the sky is also quite 
polluted - from a CIB point-of-view - by these foregrounds 
full of galactic dust. These will have to be removed from 
our maps although some residuals might remain in the 
final CIB map used for the cross-correlation. There may 
even be a significant CMB residual in this map, due to 
an imperfect separation of components. Consequently, we 
need to assess the impact of these contaminants in our study. 



As evoked in Section [2. 3[ we see that the largest contri- 
bution to the SNR comes from multipoles lower than ~ 50. 
On the other hand, the most interesting feature of these re- 
sults is that contrary to what could be intuited from Fig. [2j 
the total SNR peaks around 850 fim instead of 250 nm for 
the cross-correlation signal itself. The reason for this is ac- 
tually quite subtle : it comes from the shape of the "noise" 



term in the SNR expression in Eq. (121, as a function of i 
namely : 



[N,f{u) = {[cr{'^)f + a 



f«(z.)Cf^«)/(2^4-l) 



For all the frequencies studied here, this "noise" has roughly 
the same amplitude relatively to its corresponding "signal" : 



This is illustrated in Fig. [4] where we plotted in the left 
panel all the [Sv]^ terms with their respective maximum 
rescaled to unity. In the middle panel, we apply the same 
rescaling factor of each [S„]'^ term to the corresponding 
[A'^i,]^ term. By doing this, we can compare all frequencies 
without changing their associated signal-to-noise ratios. On 
the resulting graph, we see that a,t £ = 100 the rescaled 
noise amplitude is roughly the same, while the signal has 
the same shape at all frequencies, except for a small shift 
in £. However there is a major difference in the shape of the 
noise power spectrum from one frequency to another : its 
slope changes depending on the frequency, with the steepest 
one for Planck 850 fj.m. Therefore its amplitude goes down 
more quickly than the others as i approaches zero where 
coincidently the signal is strong, which then boosts the 
SNR at the low multipoles, and the total SNR. 

In light of these results, the optimal frequency for ISW 
detection appears to be around 353 GHz/850 fim with a 
maximum SNR reaching 7cr. However in pratice, the CIB 
extraction at this frequency might prove challenging since 
the CMB becomes dominant here, and increasingly so as 
we go down in frequency. Therefore the possible residuals in 
our extracted CIB map have to be accounted for, and other 
sources of noise as well, which is the purpose of the next 
subsection. 



3.2 More realistic SNR 

We now carry a more realistic study by including several 
possible sources of contamination : first the signal is com- 
pletely dominated on a large part of the sky by emissions 
from our own galaxy. The contamination from this fore- 
ground in the galactic plane is several orders of magnitude 
above the CIB level and prevent us from extracting the 



To account for these effects on the detectability of the 
CIB-CMB cross-correlation, we use in the present section 
a more realistic formulation of the signal-to-noise ratio, by 
adding new elements in the noise term. It therefore becomes 
at a given frequency u : 



\Cf(i. 



\Cf(y)^Nf(yy\' + [Cf B(,,)+ivCiB(^)]pCMB+ArCMB] 

(13) 

where /sky is the fraction of the sky common to the 
CMB and the CIB maps, and iVf, A^f™ and TVf^^ 
are the noise contributions respectively in the cross, CIB 
and CMB signal. Since the CMB is expected to be only 
variance-limited at the multipoles of interest, we take here 
^CMB ^ Q_ However we still have to take into account the 
CIB contamination. 

To do so, we first break the CIB noise power spectrum 
into several independant parts : 

where these four difi'erent terms represent, from left to 
right, the power spectra of the CMB residual, the galactic 
foreground residuals, the instrumental noise and finally the 
noise due to correlation between residuals and the CIB 
(which appears when autocorrelating the final CIB map). 

We quantify the CMB residual in the CIB map as a 
fraction /cmb of the total CMB map, which affects both the 
cross-correlation and CIB noise. This consequently defines 
the noise in the cross signal : 

and the following two contributions : 

rjCMB/ N f2 I \ ^, ^CMB 

Nr'°'-{u) = 2fcMB{iy)xcr{i^) 

We then define the spectrum of the foreground residuals as 
the following power law : 

so that their amplitudes are defined relatively to the real 
CIB signal through a chosen constant .4forc., which defines 
the quantity : 
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Figure 4. "Signal" terms (left-hand panel, rescaled to unity) and "noise" terms (middle panel, same rescale factor as the "signal") of 
the SNR as functions of £ (see text for details) for our chosen frequencies and instruments. The quotient of the two terms, used in the 
calculation of the SNR itself, is shown in the right-hand panel : the main difference throughout the frequencies comes from the shape of 
the "noise" term. 



i.e. the ratio between the foreground residuals and the 
CIB spectrum at the multlpole £ = 10, approximatively 
where the cross-signal is at its maximum. The slope of 
the spectrum a is fixed here for all frequencies ; previous 



analysis of infrared maps (Wright 1998 Miville-Deschenes 



et al. 20071 found it to be ~ —3 for foregrounds at high 
galactic latitudes. Finally, the instrumental noise power 
spectra A^^""*"^' at each frequency are taken from the first 



ten months of Planck data in Planck Collaboration (2011al, 



and extrapolated to the thirty months, i.e. the end of the 
fourth Planck full-sky survey. 

In this section we focus on four of the five previously 
described Planck HFI freqencies, from 217 to 857 GHz : 
we discard the fifth 143 GHz as the CMB completely 
dominates the CIB signal there. We also put aside the IRAS 
frequency here because of its weaker significance, and the 
SPIRE frequencies since the instrument is not scheduled 
to ever cover very large regions of the sky (i.e. /sky ^ 1), 



dramatically decreasing the realistic SNR (see Eq. (13l) 



At this point, we get three free parameters at each of 
the four frequencies in our SNR analysis : /sky, /cmb and 
.4forc.. The next step would be to explore this 3D parame- 
ter space at each frequency and compute the SNR at each 
point. Considering the very large number of possible combi- 
nations of parameters, it would not be practical to display 
the complete results of this exploration here. Therefore we 
first choose to fix /sky to two values of interest ; 

• ./sky ~ 0.75, which corresponds to an optimistic case 
where the only part of the sky that we discard is the galactic 
plane ; unfortunately there are other highly contaminated 
regions where the component separation techniques might 
not be able to extract the CIB. 

• /sky = 0.15, which is a low estimate of the area of the 



sky where the current data allow for an efficient CIB extrac- 
tion. The methods currently employed are based on the use 
of HI maps as a tracer of the galactic dust, though it only 
remains valid for an HI column density lower than a spe- 
cific threshold (see Planck Collaboration|[2011b for details 
on these methods). 

Concerning our other two parameters we limit ourselves to 
reasonable values, with /cmb G [0,0.1] and .4forc. G [0,10]. 

We then focus on the frequency that gave the best 
SNR results in the ideal case, namely 850 /im/353 GHz, 
and study the effect of the noise on the cross-correlation 
detectability. The results are presented in Fig. [5] which 
shows the contour levels of the SNR in the {/cmb, Aioic.) 
parameter space. The influence of the CMB is clearly 
visible at this frequency, quickly reducing the SNR as its 
residual level increases. This effect is even more pronounced 
at 1380 fim/217 GHz, where the SNR is typically twice 
as low as in the ideal case (see Table [T]), due to the fact 
that we get closer to the maximum of the SED of the 
CMB. It makes this frequency far less significant for the 
ISW detection than in the ideal case. The presence of 
instrumental noise - whose effect cannot be appreciated 
with Fig. |5] alone - becomes significant at the two lowest 
frequencies (217 and 353 GHz), again reducing their value 
in the cross-correlation. As expected the galactic foreground 
residuals also decrease the SNR, though their influence 
is roughly the same at all frequencies as they are defined 
relatively to the CIB spectrum in our analysis. Lastly, the 
biggest influence comes from the fraction of the sky through 
the /sky parameter, as the total SNR scales as \/7sky. This 
makes it a crucial requirement for future applications to 
have the largest possible coverage to minimize this effect. 
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850yLim/353GHz (f3, =0.75) 




550/2m/545GHz (f^, =0.75) 



10 



fcMB(%) 

850yLtm/555GHz (f3,y=0.15) 



10 



fcMB(%) 



Figure 5. Total signal-to-noise ratio of the CIB-CMB cross- 
correlation at 353 GHz, as a function of the CMB residuals (in 
percentage of the total CMB signal) and the foregrounds resid- 
uals (through the parameter Aioie.)- Upper panel : /gj^y = 0.75, 
the results go from less than 4 to more than 5, from the brightest 
colored area to the darkest. Lower panel : f^^^y = 0.15, the SNR 
goes from slightly less than 1.5 to more than 2, again from the 
brightest to the darkest area. 



Taking all these remarks into account and after 
some exploration of the parameter space, the optimal fre- 
quency that stands out in these more realistic scenarii is 
545 GHz/550 /im. Indeed, it is weakly influenced by in- 
strumental noise and CMB residuals and also has a higher 
"original" SNR (in the ideal case) than the other remain- 
ing frequency 857 GHz/350 /im. Our analysis at 545 GHz is 
presented in Fig. |6] 

3.3 Joint SNR 

Until now we have only considered a detection at a single 
CIB frequency and its associated significance. In practice, 
we will have several cross spectra at different frequencies, 
e.g. in the case of Planck where we will be able to extract 
the CIB at four different frequencies on a large fraction of 
the sky. This allows us to increase the total signal-to-noise 
ratio of the ISW detection by combining the constraints 
from all available frequencies, though this will be limited by 
the possible intrinsic correlations between the CIB maps at 
different frequencies. Indeed, such correlations imply some 




fcMB(%) 

550/xm/545GHz (1^0.15) 




Figure 6. Total signal-to-noise ratio of the CIB-CMB cross- 
correlation at 545 GHz as a function of the CMB residuals and 
the foregrounds residuals. Upper panel : /gjjy = 0.75, the results 
go from slightly less than 4 to more than 5, from the brightest to 
the darkest area. Lower panel : /gj^y = 0.15, the SNR goes from 
less than 2 to slightly more than 2.5. 



redundancy in the measured information, and therefore 
lessen the gain in the total significance of the combined 
detection. 

We can expand the previous SNR formalism to ex- 
press the theoretical joint significance of a set of n cross- 
correlations (i.e. CIB at n frequencies, each correlated to 
the same CMB) : 



(14) 



with X {X'^) being the column (row) vector of all the cross- 
correlations : 



x^ = (x^(i^i) ••• x^K)) 



where X'^{ui) contains the cross-spectrum at the frequency 
from ^ = 2 to 100 : 

X^i^i) = (CeL^iui) ■ ■ ■ CeUooii^i)) 

The block matrix A4 is the covariance matrix, containing 
n X n blocks. Each one of them represents the covariance of 
two cross-spectra at different CIB frequencies, depending on 
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the position of the block. At the ith line and jth column, 
the block A4''-' is written as : 



\ 







The diagonality of M^-' comes from the assumption that the 
different multipoles are uncorrelated. In the noiseless case 
discussed in Section [3. 1[ the elements of each block can be 
expressed as follows : 

M;^' =Covar(C?(z.i),C7(j^j)) 
21+1 

We can see here the dependence on the aforementionned 
possible correlation between the CIB at frequency Ui and 
the CIB at frequency Vj, through the cross-spectrum 
Cf'^^^{vi,Vj). To perform a more advanced analysis, we 
can easily modify this expression to account for the possible 
sources of noise discussed in the previous section. 

Once again, the large number of possible combinations 
of noise parameters makes it unpractical to present a 
complete study of the joint correlation. Instead we focus 
on a few particular cases, motivated by what we found 
in Section |3.2| A summary of our results on single and 
joint correlation is presented in Table [l] : we first go 
back to the ideal case to quantify the impact of the joint 
detection. We found a relatively small gain, as it increases 
the total SNR by a mere ~ 0.15 compared to the maximum 
significance of a single detection. This can be attributed 
to the high correlations between the CIB at its different 
observed frequencies, which limit the usefulness of the joint 
cross-correlation. 

Considering now more realistic situations, with the 
presence of instrumental noise, we once again choose to fix 
some of the parameters mentionned in Section |3.2| with 
/sky = 0.75 and /sky = 0.15. A reasonable confidence 
in compenent separation techniques allows us to hope for 
small enough residuals, so that we choose /cmb = 0.01 and 
^forc. = 0.01. In these cases, the joint correlation has once 
again a limited interest (respectively a ~ 0.15 and ~ 0.07 
gain for /sky = 0.75 and 0.15) due to the correlations in both 
the CIB signals but also in the astrophysical noise contribu- 
tions - CMB and dust - between frequencies. 



4 CONCLUSIONS 

The topic of this paper is an investigation of the cross- 
correlation between the cosmological infrared and microwave 
backgrounds, and a study of its detectability under various 
observational situations. A non-zero correlation is expected 
to exist between the two backgrounds and their anisotropies 
through the ISW effect, caused by the time-evolving grav- 
itational potentials that underlie the large-scale structures 
which are the sources of the CIB and of its anisotropies. 
Describing the CIB anisotropies as linearly biased tracers of 



Frequency (GHz) 


857 


545 


353 


217 


Wavelength (fim) 


350 


550 


850 


1380 


Perfect Single SNR 


6.26 


6.83 


6.98 


6.95 


Joint SNR 


7.12 


Realistic Single SNR 1 

(/sky = 0.75, /cmb = 0.01, 

Afore. = 0.01) 


5.36 


5.73 


5.39 


3.56 


Joint SNR 


5.88 


Realistic Single SNR 2 

(/sky = 0.15, /cmb = 0.01, 

Aiorc. = 0.01) 


2.40 


2.56 


2.41 


1.59 


Joint SNR 


2.63 



Table 1. Total signal-to- noise ratio of the CIB-CMB cross- 
correlation for four of the CIB frequencies of PZancfc-HFI. The 
results are given for each frequency and for the joint cross- 
correlation, first for the ideal case discussed in Section |3.1| and 
then for two more realistic cases. 



the matter field fiuctuations, we calculated the theoretical 
angular power spectrum of the CMB-CIB cross-correlation 
at several frequencies and for different instruments, taking 
into account their actual bandpasses. As is well known for 
CMB/galaxies cross-correlations, the signal peaks at low 
multipoles and quickly vanishes at higher I. The linear 
bias introduced by our formalism was then obtained by 
confronting our predicted linear CIB power spectra with the 
data coming from the Planck mission. These observed CIB 



anisotropies were fitted in Planck Collaboration (2011b I by 



an HOD model, to which we compared our own spectra at 
the low multipoles in order to get the desired bias at each 
different frequency. 

Using an advanced SNR analysis which included the 
main sources of noise both instrumental and astrophysical, 
and all their possible correlations, we pointed out the most 
promising frequency in the ideal case of noiseless full-sky 
maps (850 /j,m/353 GHz) with an expected significance as 
high as ~ 7(7 for the cross-correlation signal. The same 
frequency turned out to be less optimal with more realistic 
assumptions about sky coverage and possible sources of 
noise (here CMB, dust residuals and instrumental noise). In 
this case, higher frequencies such as PZancfc-HFI's 545 and 
857 GHz are favored, with an expected significance ranging 
from 2.4 to 5.7 depending on the frequency, the levels of 
noise and the fraction of the sky available for analysis. We 
also found that a joint cross-correlation using all available 
frequencies is of minor interest, due to the high correlations 
between CIB anisotropies at the different frequencies. 
Nevertheless, our best results for /sky = 0.75 are higher 
than the significances of all current CMB-galaxies cross- 
correlation, with (T > 5, although a less optimistic estimate 
for the sky coverage quickly reduces our signal-to-noise 
ratios. This stresses once again the requirement of good 
component separation techniques and foreground removals 
for future applications, in order to have the largest fraction 
of common clean sky /sky possible. 

The results of this work will be valuable in the forth- 
coming years of analysis and exploitation of the Planck 
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data. The formalism we developped provides us with an 
accurate forecast of the expected results of the CIB-CMB 
cross-correlation and allows us to constrain the requirements 
for a significant ISW detection. Regarding the use of the 
CIB itself, it presents some advantages over classical ISW 
studies : the underlying structures observed through the 
CIB span a large integrated range of redshifts and cover the 
whole sky whereas the usual galaxy surveys used often have 
a limited depth and width in redshift or a small sky cover- 
age - some of the main limiting factors in the ISW detection. 

In current studies, the CIB is always considered in its 
integrated form at a given frequency, meaning that in this 
observed CIB, contributions from many redshift ranges are 
mixed together. An interesting further step would be to 
use the multiple observed frequencies to reconstruct the 
contributions from different redshift bands, in order to ob- 
tain several decorrelated CIB maps corresponding to these 
redshift slices. The resulting independent CIB maps could 
then be individually correlated with the CMB. Combining 
the independent detections could increase even more the 
total SNR of the ISW detection, as it allows to get rid 
of the correlation terms between CIB maps. Fiirthermore, 
each of these maps will then help tracing the Dark Energy 
at a different time. Our preliminary calculations from 
predicted power spectra indicate encouraging enhancements 
in the signal-to-noise ratio, although the details of the CIB 
decorrelation need further investigation and optimization, 
and will be presented in a future work. 

Finally, let us note that as a background the CIB is 
likely to be lensed by large-scale structures in the local Uni- 
verse : a dedicated study of the effects of lensing in a future 
work will be able to determine if the lensing could lead to a 
possible gain in the signal-to-noise ratio of the ISW effect, 
or should be considered as a possible source of bias in the 
DE detection. 
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